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Calculations are presented for electronic Raman scattering and infrared conductivity in a d x 2_ y 2 
superconductor including the effects of elastic scattering via anisotropic impurities and inelastic 
spin-fluctuation scattering. A consistent description of experiments on optimally doped Bi-22f 2 is 
made possible by considering the effects of correlations on both inelastic and elastic scattering. 
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Impurity effects in the cuprates have played a major 
role in clarifying the nature of unconventional supercon- 
ductivity. Via the change in the density of states (DOS) 
near the Fermi level, crossover behavior of low tempera- 
ture transport and thermodynamics quantities could be 
explored as a testing ground for unconventional pairing 
0. Examples include the T to T 2 crossover in the low 
temperature magnetic penetration depth A(T), T 3 to T 
crossover in the NMR rate, and the w 3 to u> crossover in 
the low frequency B\ g Raman response. Simple theories 
based on nearly-resonant scattering in the impurity T— 
matrix approximation ("dirty" d-wave) have accounted 
for these behaviors seen in materials of various quality 
and impurity dopings. 

However, several inconsistencies arise when attempt- 
ing to put together a complete picture. First, the mag- 
nitude of the impurity scattering needed to fit A(T) and 
the extrapolated T = resistivity are generally much 
smaller than those needed to fit the measured frequency 
dependence of the infrared conductivity (IR) and the 
Raman response. In particular, the impurity scattering 
rate 1/Tj mp = 2T needed to fit the IR on high quality 
Y-123 was 100 times larger than that needed to fit mi- 
crowave measurements 0], while values of T/Ao rang- 
ing up to 0.5 were needed to interpret the Raman data 
taken on Bi-2201 |§, and some samples of Tl-2201 @] 
and Hg-1223 Q. Second, the "universal" dc conductivity 
a(to — > 0, T = 0) = (e 2 /2irTi)t;/ a || has not been convinc- 
ingly observed and the low temperature variation of the 
conductivity Sa(T) changes slower than T 2 [Q. Third, 
the transition temperature T c is only moderately reduced 
by planar impurities compared to Abrikosov-Gorkov the- 
ory H (by about a factor of two to three). Lastly, the 
residual resistivity obtained for Zn doped Y-123 suggests 
that a strong contribution is present in d— wave scatter- 
ing channels [||. 

These inconsistencies reveal that the usual treatment 
of point-like s— wave impurities in a T— matrix approach 
may be too naive and neglects electronic correlations. 
Experimental evidence from transport measurements in 
Zn doped YBCO have revealed that the a single impu- 
rity embedded in the Cu02 plane disturbs the local en- 



vironment and yields an effective scattering cross sec- 
tion diameter of Zn 2+ of 4.2A |l(J and is thus extended. 
Theoretical studies have shown that static vacancies in a 
Heisenberg antiferromagnet enhance the staggered mag- 
netic moment within a few lattice spacings around the va- 
cancy |ll[] , while studies of models with only short range 
antifcrromagnetic (AF) correlations have shown that cor- 
relations dynamically generate finite-range impurity po- 
tentials from single site impurities [Q. 

In this paper we explore the role extended impurities 
have on transport properties in an attempt to resolve 
the abovementioned discrepancies. In particular, we re- 
examine the IR and the electronic Raman response of 
Bi-2212 including the effects of electronic correlations on 
both inelastic and elastic scattering potentials, and com- 
pare our results with the extrapolated T = normal state 
resistivity and the crossover temperature for A(T). We 
find that a consistent picture emerges when we include 
the extended range of impurity scattering as well as AF 
spin fluctuations. 

Currently, knowledge of a T— matrix formulation for 
disorder in correlated systems is limited and generally ap- 
proximate methods have been used fll3|] . Formally, one 
would need to include not only bare impurity and inter- 
action self energies responsible for elastic and inelastic 
scattering, respectively, but also one must include terms 
in which mix elastic and inelastic potentials. In this 
way, a purely static bare impurity interaction changes 
in nature (can become dynamic) due to the inclusion of 
many-body effects. Ziegler et al. in Ref. (T^| discussed 
how a point-like s— wave bare impurity potential can be- 
come extended due to the background of the correlated 
host. One can then proceed to calculate the T— matrix 
by using a renormalized Hamiltonian which describes the 
impurity potential in the correlated host. 

Our starting point is the model Hamiltonian consid- 
ered in Ref. fll4| ] , which represents the effects of impuri- 
ties in a metal with short range AF order on a square 
lattice: 
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The parameters Vq and V\ denote the on-site and ex- 
tended impurity potentials, and tj denotes the effect of 
impurities on the electron hopping to the impurity site. 
Focusing on a two parameter model using the specific re- 
lation V\ — a 2 Vb/4 and tj — aVb/4, an analytic solution 
for the single impurity T-matrix was obtained, a is the 
control parameter which distinguishes between point- like 
(a = 0) and extended (a ^ 0) impurity potentials. 

The algebraic solution for the impurity T-matrix is 
presented in Ref . jL4j , where the reader is referred to for 
details. In essence, the extended structure of the impu- 
rity potential requires a 4 x 4 matrix formulation with 
respect to s, p, and d scattering channels expressed in 
terms of the k dependent basis functions of the square 
lattice. The impurity averaged self energy is determined 
via £* mp (k, iuj) — riitk^iu)), with rij the impurity con- 
centration. The self energy can be expanded in Pauli 
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and S™ p (k, iui) = 0. Here so,3, <io,3,Po,3, 8, (3 and a* 
are constants [|15| given by Eqs. (8-11) of Ref. Q . 
The k-dependent basis functions for a square lattice 
are 7k, 7k = [cos(k x a) ± cos(k y a)}/2, and 7^, j^ 2 = 
[sin(fc^a) ± s'm(k y a)]/2. In the limit of a — > 0, the T- 
matrix becomes momentum independent and recovers 
the well known results . 

In the self consistent T— matrix approximation, the 
non-interacting Green's function is used to calculate 
the self energy and then the new Green's function is 
put back into the self energy calculation. The pro- 
cess is iterated until convergence is realized, which oc- 
curs typically after only a few iterations. We assume 
a strong on-site impurity interaction Vq = 8t, take 
£k = —2t{cos(k x a) + cos(k y a)] + At' cos(k x a) cos(k y a) — p, 
Ak(T) = A (T)[cos(fca;a) — cos(k y a)]/2, a weak coupling 
form for A (T) and choose A (T = 0) = 4T C = OAt. 
Here and throughout, lattice sizes of 32 x 32 up to 
256 x 256 and typically over 1000 frequencies were used. 
Our results showed little size effects above the 128 x 128 
mesh. The chemical potential fi was adjusted so that the 
filling (n) = 0.825. 

In the absence of vertex corrections, the homogeneous 
Raman response and the real part of the conductivity are 
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FIG. 1. (a) Bi g (upper set) and Bi g (lower) Raman re- 
sponse for T = 0.2T C , (7 = and 1% impurities. Here and in 
(b) the solid (dotted, dashed, long-dashed) line corresponds 
to a = 0(0.5,0.75, 1). (b) Log-log plot of the low frequency 
Bi g Raman response determining the crossover frequency ui* . 
The legend gives the values of lu* /t. 
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where / is the Fermi function and the current vertex 
ik = e de\s_/dk x = 2teasm(k x a){l — 2t' /tcos(k y a)]. The 
Raman response can be classified according to elements 
group: 
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If the light scattering is non-resonant, the frequency de- 
pendence of the momentum independent prefactors b, b' 
and a can be safely neglected, and we can adjust these 
prefactors to account for overall intensity. It can be seen 
from the k— dependence of the vertices that the B\ g re- 
sponse probes qp dynamics around the Brillouin Zone 
(BZ) axes, B2 g probes the diagonals, and A\ g is more 
of an average around the BZ and involves pure density 
fluctuations and backflow jl7| . Higher order terms of in- 
creasingly more anisotropic basis functions could be con- 
sidered but do not lead to appreciable differences except 
for the Ai g response (see Ref. jl^] for details). From here 
on we only consider the B\ g and Bi g channels. 

Results for the B\ g and B2 g response are plotted in 



Fig. (|l|a) for different values of a. As seen in our previous 
studies of the DOS (lij, turning on a even slightly leads 
to an effective increase in the strength of the impurity po- 
tential. The increased impurity resonance at the Fermi 
level is manifest in the increasingly smeared spectra as 
well as the value of u*, defined as the frequency where 
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FIG. 2. The impurity scattering rate for momenta along 
the Fermi surface for a = 0, 0.24, 0.48, and 0.72 (solid, dotted, 
dashed, long-dashed lines, respectively). 
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FIG. 3. The calculated IR at T=0.2T C , (7 = 0, and 
2% impurities for a = 0,0.25,0.5,0.75 and 1 (solid, 
dashed, long-dashed, dash-dotted lines, and thin solid 
lines, respectively) . 



the cubic frequency dependence becomes sub-dominant 
p9[ , as shown in Fig. (0b). Since a = 0.5 corresponds 
to a nearest neighbor impurity potential V± = Vb/16, the 
presence of even a 6% nearest neighbor potential leads to 
a 23% increase in lu*. Therefore a much smaller concen- 
tration of extended impurities is needed to have the same 
effect as isotropic impurity scattering. The B2 g channel 
is only slightly affected by a. 

This is in accord with the impurity scattering being 
largest near the BZ axes in our model and only minorly 
dependent on a near the BZ diagonals. This is seen in 
Fig. (JsJ) , which plots the normal state T — impurity 
scattering rate around the Fermi surface for different a. 
Increasing a from zero first decreases then increases the 
scattering rate near the BZ axes while only mildly affect- 
ing the rate along the BZ diagonals. Therefore we would 
expect the B\ g channel to be more sensitive than Bi g to 
the growth of scattering for increased a. 

As an important consequence, the IR (which has sim- 
ilar weighting to i?2g) does not pick up regions where 
the scattering changes rapidly and is governed by small 
scattering along the BZ diagonals. The IR in the super- 
conducting state is shown in Fig. (||) for different values 
of a at a fixed impurity concentration. The IR does not 
appreciably change for a up to 0.5. For larger a spec- 
tral weight is shifted away from low frequencies as the 
scattering rate increases over the entire BZ |po| , and a 
stronger shoulder near A max develops. 

We now consider inelastic scattering in order to address 
the intensity observed at high frequencies. The origin of 
the flat background in IR and Raman at high frequen- 
cies has attracted a great deal of attention involving both 
Fermi liquid and non- Fermi liquid based theories. A cor- 
rect approach would account not only for the featureless 
continuum but would be able to describe the polarization 
dependence of the Raman cross section and the differ- 
ences between IR and Raman. Raman has shown strong 
evidence of two-magnon features in the insulating as well 



as superconducting state, implying the role of spin fluc- 
tuations as a source for inelastic scattering |2l|. The 
way in which spin fluctuations are included in calculat- 
ing dynamic quantities has also attracted a large amount 
of attention. The main problem has been the degree in 
which strong local electron correlations arc included and 
represented. Effective Hamiltonian |22] approaches based 
on Bragg-like scattering of quasiparticles for momentum 
Q = (-7T, 7r) have been successfully employed to calculate 
various dynamic correlation functions, including IR ]23j ] 
and Raman j24|. However, to a certain extent these ap- 
proaches have their limitation as they do not adequately 
capture the strong multi-magnon scattering process re- 
quired to approach the insulating phase from the metallic 
side [ p5[ . 

Since here we are interested in only the low frequency 
behavior of the IR and Raman response, the details of the 
dynamic scattering are not crucial in that they only affect 
the response functions at larger frequencies. Therefore 
we take a simple route and include spin fluctuations in 
RPA, 
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where U is a phenomenological parameter [we choose U = 
Xoi^i^) is the non-interacting spin susceptibility, 
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FIG. 4. The IR Conductivity for T = 0.2, 0.8, 0.9, and T c 
(solid, dashed, long-dashed, dash-dotted lines, respectively) 
for 2% impurities and a = 0.5 and U — 2t. The thin solid 
line is for U — 0, T = 0.2T C for comparison. 
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The convolution of momentum sums in Eq. ^ is solved 
numerically via Fast Fourier Transform, where we keep 
the full k-dependence and the real and imaginary parts of 
the self energy. We found that neglecting the real parts 
of the self energy and/or restricting momentum sums 
around the FS leads to a substantially smaller conduc- 
tivity and misses a renormalization of the conductivity 
peak to frequencies slightly away from 4A in the super- 
conducting state Combining both Y, lmp,u , the results 
for the IR and the Raman response for the B\ g and Bi g 
channels are summarized in Figs. (Q) and (0). For both 
quantities the spin fluctuations yield a flat continuum at 
high frequencies in common with experiments and vari- 
ous different theories which yield a linear frequency de- 
pendence of the imaginary part of the self energy. The 
temperature dependence there is minimal and all the 
spectra converge to similar values by roughly f2 ~ 2t. 
As the temperature is lowered, the low frequency IR falls 
in magnitude and develops a shoulder at ~ A max and a 
weak peak at ~ 4A maa; , while the spectral weight in the 
Big and i?2g channels reorganize from low frequencies to 
higher frequencies at ~ 2A max and 0.65A max , respec- 
tively. Strong peaks associated with pair breaking and 
the van Hove structure (~ t) appear in the B± g channel 
and become less pronounced as the temperature is in- 
creased due to the growth of spin fluctuation scattering. 
We note that at present there is no experimental indica- 
tion of a peak in the B\ g which could be associated with a 
van Hove feature J3| . In our calculations further smearing 
of the van Hove peak is expected if dispersion is added in 
the c— direction or if stronger interactions are used which 
produce larger inelastic scattering at large frequencies. 



FIG. 5. Raman response for the B lg (a) and B29 (b) chan- 
nel for U — 2t, Vo — St, 2% impurities, and a = 0.5 for 
T/Tc = 0.2, 0.6, 0.8, 0.9, 1 plotted as the thick solid, dotted, 
dashed, long-dashed, and dotted-dashed lines. The thin solid 
line shows the response at T = T c for U = for comparison. 



Moreover, a more correct multi-magnon approach would 
also wash out structure at higher frequencies. 

Fits to the Raman spectra measured in the supercon- 
ducting state at T = 0.5T C of an as-grown sample of 
Bi-2212 are presented in Fig. ||, while fits to the IR in 
both the normal and superconducting state on a similar 
sample with a slightly higher T c by taken by N. L. Wang 
et al. |^| are shown in Fig. [?]. Here we have taken the 
parameters used in Fig. have adjusted the prefactors 
bjb 1 = 1.46 to account for the relative Raman intensities. 
For the IR, we use the c— axis lattice spacing of 30A for 
Bi-2212 containing 2 Cu02 bilayers to convert the 2D IR 
to 3D. We find that theory underestimates the IR scale 
by only a factor of 1.5 (the fits in Fig. |t] are scaled by 
this factor). We used slightly different values of t for Ra- 
man (81meV) than IR (69 meV), but equally good fits 
are obtained if we used different values for A. The results 
agree exceptionally well with the measured spectra espe- 
cially at low frequencies where the effects of impurities 
are dominant. The agreement lessen to to only a quali- 
tative level for > 1000 cm -1 due to the small degree 
of spin fluctuations included (in RPA) and points to an 
inadequate description of the normal state. 

Finally, we consider how the fitting parameters com- 
pare to transport data. Assuming a Drude model for 
p(T = 0)(~ lO/iO-cm) and a plasma frequency of 1.2 
cV [|?J implies a scattering rate l/T imp = 15cm" 1 , while 
the penetration depth T to T 2 crossover [ p8| measured in 
the same sample |2j| as the Raman data gives 12 cm -1 . 
Previous Raman fits using isotropic impurity scattering 
only and a Fermi surface restricted approach required 
l/Timp = 2r = 72 cm -1 |1^| . Our calculation for a = 0.5 
yields 1/t$£ = 23 and 20 cm" 1 for the FS averaged 
impurity scattering rate for the t = 81 and 69 meV, re- 
spectively. Given the uncertainty in estimating p(T = 0) 
and ujpi, this in favorable agreement with existing mea- 
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FIG. 6. Fits to the Bi g (a) and B 2a (b) low temperature 
spectra on Bi-2212 (T C =86K) taken by R. Hackl et al. in ||. 
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FIG. 7. Fits to the normal and superconducting IR on 
Bi-2212 (T C =93K) taken by N. L. Wang et al. [§|. 

surements. 

In summary we have shown how the inclusion of elec- 
tron correlations in both inelastic and elastic scattering 
potentials lead to a consistent description of the channel 
dependent Raman and IR response and lead to better fits 
than previously achieved. The overall intensity of the IR 
can be accounted for and the effect of extended impuri- 
ties on the low frequency behavior of the Raman response 
can resolve the discrepancy between large impurity scat- 
tering rates needed previously for IR and Raman fits and 
the small rates needed for transport. 
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